New method of precise measurement of positronium hyperfine splitting 



A. Ishida 3 -*, G. Akimoto a , Y. Sasaki 3 , T. Suehara\ T. Namba\ S. Asai a , T. Kobayashi b , H. Saito c , M. Yoshida d , K. Tanaka d , and 

A. Yamamoto d 

"Department of Physics, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
b International Center for Elementary Particle Physics (ICEPP), The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
' Department of General Systems Studies, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro-ku, Tokyo 153-8902, Japan 
" ! High Energy Accelerator Research Organization (KEK), 1-1 Oho, Tsukuba, Ibaraki 305-0801, Japan 



Abstract 



o 

(N 

5— ( 
Oh 

< 

o 



The ground state hyperfine splitting of positronium, Ahfs, is sensitive to high order corrections of QED. A new calculation up 
to OCc 3 In a) has revealed a 3.9<x discrepancy between the QED prediction and the experimental results. This discrepancy might 
either be due to systematic problems in the previous experiments or to contributions beyond the Standard Model. We propose 
an experiment to measure Ahfs employing new methods designed to remedy the systematic errors which may have affected the 
previous experiments. Our experiment will provide an independent check of the discrepancy. The prototype run has been finished 
and a result of Ahfs = 203.3804 + 0.008 4 GHz (41 ppm) has been obtained. A measurement with a precision of O(ppm) is expected 
within a few years. 
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1. Introduction 

Positronium (Ps), a bound state of an electron and a positron, 
is a purely leptonic system which allows for very sensitive 
tests of Quantum ElectroDynamics (QED). The precise mea- 
surement of the hyperfine splitting between orthopositronium 
(o-Ps, l 3 5i) and parapositronium (p-Ps, l'So) (Ps-HFS) pro- 
vides a good test of bound state QED. Ps-HFS is expected to 
be relatively large (for example compared to hydrogen HFS) 
due to a relatively large spin-spin interaction, and also due to 
the contribution from vacuum oscillation (o-Ps — > y* — > o-Ps). 
The contribution from vacuum oscillation is sensitive to new 
physics beyond the Standard Model. 

Figure Q] shows the measured and theoretical values of Ps- 
HFS. The combined value from the results of the previous 2 
experiments is A^ = 203.388 65(67) GHz (3.3 ppm) flU. 
Recent developments in NonRelativistic QED (NRQED) have 
added 0(or 3 In a) corrections to the theoretical prediction which 
now stands at A* FS = 203.391 69(41) GHz(2.0ppm) [3]. The 
discrepancy of 3.04(79) MHz (15 ppm, 3.9cr) between A^? s and 
A* FS might either be due to the common systematic uncertain- 
ties in the previous experiments or to new physics beyond the 
Standard Model. 

There are two possible common systematic uncertainties 
in the previous experiments. One is the unthermalized o-Ps 
contribution which results in an underestimation of the mate- 
rial effect. This effect has already been shown to be signifi- 
cant 10, El 0] in the o-Ps lifetime puzzle. The other is the uncer- 
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Figure 1: Measured and theoretical values of Ps-HFS. 



tainty in the magnetic field uniformity which was cited as the 
most significant systematic error by previous experimenters. 

2. Theory of Experiment 

2.1. Measurement using Zeeman effect 

The energy levels of the ground state of Ps are shown as a 
function of static magnetic field in Figure [2] Due to techni- 
cal difficulties in directly stimulating Ahfs, we make an indi- 
rect measurement by stimulating the transition A m ; x . This is the 
same approach as previous experiments. The relationship be- 
tween Ahfs and A m ; x is approximately given by the Breit-Rabi 
equation 



A mix * iA HFS (Vl+4x 2 -l) 
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From the time-dependent Schrodinger equation, the 4x4 den- 
sity matrix p(t) is given by 



ihp = Hp- pH f 



(4) 



where the i,j-element of pit) is defined as pij(f) = 
(iffi\i/r(t))(iff(t)\iffj). If we take the initial state to be unpolarized, 
p(0) = diag(i,i,i,i). 

The 2y-ray annihilation probability S 2 r and 3y-ray annihila- 
tion probability Sj, y , between t = to and t = f], are obtained 
by 



Siy = r, 
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Figure 2: Zeeman energy levels of Ps in its ground state. The arrows f, I means 
the spin up and down of electron, and the arrows ft, li means the spin up and 
down of positron. 



in which x = g'psB/hAuFs, g' = g (l - ^£* 2 ) is the g factor for 
a positron (electron) in Ps 07[], ps is the Bohr magneton, B is 
the static magnetic field, and h is the Plank constant. 

In a static magnetic field, the \S,m z ) = |0, 0), where S is 
the total spin of Ps and m z is the magnetic quantum number of 
Ps along with z-axis (direction of static magnetic field), state 
mixes with the |1,0) state hence the |+) state annihilates into 
2 y-rays with a lifetime of about 8 ns (with our experimental 
conditions). The |1,+1) states annihilate into 3 y-rays with a 
lifetime of about 140 ns. When a microwave field with a fre- 
quency of Amix is applied, transitions between the |+) state and 
the |1,+1) states are induced so that the 2 y-ray annihilation 
rate increases and the 3 y-ray annihilation rate decreases. This 
change of annihilation rates is our experimental signal. 

Our experimental resonance line shape is obtained using den- 
sity matrix. We use the basis for four spin eigenstates of Ps as 
0Ao,<Ai,iA 2 ,<A3) = (|0,0>,|1,0),|1,1),|1,-1». We apply a mag- 
netic field 



B(t) = Be z + B e x cos (cot) , 



(2) 



where e z , e x are the unit vectors for z, x direction respectively, 
Bq is magnetic field strength of microwaves, and a> is the fre- 
quency of microwaves. Then the Hamiltonian H becomes 



H 



-I - Ly 

2 2/s 

A' 

-y 

y 



7t 



o 
o 



-y 
o 

5-ir. 
o 



y 
o 
o 



>7t J 



/*A HFS , (3) 



'p-Ps 

2ttAhfs : 
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where y = C y 8 h ^ n> ° cos (u)t), C y is a constant, y s 

2^^ » Tp-Ps is the decay rate of p-Ps, and T _p s is that of o-Ps. 
The most recent and precise experimental values are r„p s = 
7.040 1(1) ps- 1 H andr p _ Ps = 7.990 9(17) ns" 1 H. 



2.2. Ps Thermalization Effect 

Forming Ps needs material to provide electrons, but material 
around Ps makes electric field and changes Ahfs- This effect 
is called the Stark effect. The material effect must be properly 
considered to evaluate Ahfs in vacuum. In the previous exper- 
iments, the material effect on Ahfs was considered to be pro- 
portional to the material (gas) density. Ahfs was measured at 
various gas density, and they were extrapolated linearly to zero 
density |jl 0]. But this extrapolation method can make large 
systematic error enough to account for the discrepancy between 
Ahfs ™d A£ FS . 

Formed Ps has the initial energy of O(eV). Positronium loses 
its energy when it collides with material, and finally its energy 
becomes to room temperature (~ 1/30 eV). This process is 
called thermalization. The material effect between t - to and 
t = t\ is proportional to 



f 

Jto 



/(f) (poo(t)(B * 0) - poo(t)(B = 0)) dt , 



(7) 



where /(f) ~ ncrv(t) is the collision rate of Ps with material, n 
is the number density of material, cr is the typical cross section 
of collision, and v(t) is the mean velocity of Ps. In the previous 
experiments, timing information was not measured so that to 
and f i = oo. The effect is proportional to the material density 
if the thermalization occurs much faster than the o-Ps lifetime, 
but the thermalization time scale becomes large especially at 
low material density, which makes nonlinear effect on Ahfs- 

The nonlinear effect can be estimated using the thermaliza- 
tion model [9], 



dEw(t) 
dt 



■V2w Ps £ av (f)(£ av (f) - ^kTj ^ 
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where £ av (f) is the average Ps energy, mp s is the Ps mass, M is 
the mass of the gas molecule, T is the temperature of the gas, k 

is the Boltzmann constant, and | -^^^7 means the collision 
effect. The solution of this equation is II 1 011 
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Figure 3: (Color online) Thermalization effect on Ahfs- The circles and error 
bars are the data of Ref. [1], the solid line is the linear fit, the dashed line is 
the best fit including Ps thermalization effect, and the red band is the 0(a 3 In a) 
QED prediction in vacuum. 



where b = A = ^j=, and E , 

E- dv {0) is the initial energy of Ps. 

Figure [3] shows the best fit result of the previous experimen- 
tal data from Ref. |Q]] considering nonlinear material effect us- 
ing Equations {7]i and (|9]l with parameters Eq = 2.07 eV and 
cr m = 13.0 A 2 11 111 , compared with linear fitting. The nonlin- 
ear effect can be clearly seen at low density, and it can explain 
the discrepancy between A™p S and A* FS (the best fit value in- 
cluding nonlinear effect is Ahfs = 203.392 80(95) GHz). But 
cr m = 37 + 10 A 2 from Ref. |0], which is not consistent with the 
value of Ref. Ill ill , affects only about 5 ppm and cannot explain 
the discrepancy. Therefore, new independent measurement of 
Ps thermalization is needed. 



2.3. Our new methods 

Our new methods will significantly reduce the systematic er- 
rors present in previous experiments; the Ps thermalization ef- 
fect and the non-uniformity of the magnetic field. The main 
improvements in our experiment are the large bore supercon- 
ducting magnet, /^-tagging system, and high performance y-ray 
detectors. Details are discussed in the following sections. 

3. Prototype Run 

The prototype run of the measurement with new methods has 
been performed. 

3.1. Experimental Setup 

A schematic diagram of the experimental setup of the proto- 
type run is shown in Figure [4] 

3.1.1. Large bore superconducting magnet 

A large bore superconducting magnet is used to produce the 
magnetic field B ~ 0.866 T which induces the Zeeman split- 
ting. The bore diameter of the magnet is 800 mm, and its length 




Figure 4: Schematic diagram of the experimental setup of the prototype run 
(top view in magnet). 



is 2 m. The large bore diameter means that there is good uni- 
formity in the magnetic field in the region where Ps is formed. 
Furthermore, the magnet is operated in persistent current mode, 
making the stability of the magnetic field better than + 1 ppm. 

The magnetic field distribution was measured using a pro- 
ton NMR magnetometer (ECHO Electronics EFM-150HM- 
AX) with a clock synthesizer (NF Corporation CK1615 PA- 
001-0312). Calibration uncertainty of the system is less than 
1 ppm. Non-uniformity of the magnetic field in a whole vol- 
ume of RF cavity is 23.1 ppm (RMS) without any compen- 
sation coils. The non-uniformity decreases to 10.4ppm when 
weighted by RF power and positron stop position distribution. 

3.1.2. B-tagging system and timing information 

The positron source is 19/vCi (700 kBq) of 22 Na (Eckert & 
Ziegler POSN-22). A plastic scintillator (NE102A) 10 mm in 
diameter and 0.2 mm thick is used to tag positrons emitted from 
the 22 Na. The scintillation light is detected by fine mesh pho- 
tomultiplier tubes (PMT: HAMAMATSU H6614-70MOD) and 
provides a start signal which corresponds to the time of Ps for- 
mation. The timing resolution is 1.0 ns (l<x). The positron then 
enters the microwave cavity, forming Ps in the gas contained 
therein. 

Ps decays into photons that are detected with LaBr3 (Ce) 
scintillators (Saint-Gobain BrilLanCe™ 380). Accumulating 
measurements of the times of positron emission and y-detection 
results in decay curves of Ps as shown in Figure [5] The timing 
information is used to improve the accuracy of the measurement 
of Ahfs as follows: 

1. Imposing a time cut means that we can select well ther- 
malized Ps, reducing the unthermalized o-Ps contribution. 
It should also be possible to precisely measure the contri- 
butions of unthermalized o-Ps, and of material effects (we 
plan to make such measurements in future runs). 

2. A time cut also allows us to avoid the prompt peak (con- 
tributions of simple annihilation and of fast p-Ps decay), 
which greatly increases the S/N of the measurement (by 
about a factor of 20). 
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3.1.3. High performance y-ray detectors 

Six y-ray detectors are located around the microwave cavity 
to detect the 511 keV annihilation y-rays. LaBr3 scintillators, 
1 .5 inches in diameter and 2 inches long are used. The scintilla- 
tion light is detected by fine mesh PMT through the UVT light 
guide. Without light guide, LaBr3 scintillators have good en- 
ergy resolution (4% FWHM at 5 1 1 keV) and timing resolution 
(0.2 ns FWHM at 5 1 1 keV), and have a short decay constant 
(16 ns). The good energy resolution and the high counting rate 
of LaBi - 3 results in very good overall performance for measur- 
ing Zeeman transitions. In particular the good energy resolu- 
tion allows us to efficiently separate 2y events from 3y events, 
negating the need to use a back-to-back geometry to select 2y 
events, thus greatly increasing the acceptance of our setup. 

This y-ray detector system greatly reduces the statistical er- 
ror in the measurement. 
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Figure 5: (Color online) Decay curves of Ps. The red line is on-resonance 
(0.865 816 9 T) RFON, and the black line is RFOFF. The timing window of 
35-155 ns is also indicated. The decay rate of Ps increases with RF because of 
the Zeeman transition. 



3.1.4. RF system 

Microwaves are produced by a local oscillator signal gener- 
ator (ROHDE & SCHWARZ SMV 03) and amplified to 500 W 
with a GaN amplifier (R&K A2856BW200-5057-R). 

The microwave cavity is made with oxygen-free copper; 
the inside of the cavity is a cylinder 128 mm in diameter and 
100 mm long. The side wall of the cavity is only 2 mm thick 
in order to allow the y-rays to efficiently escape. The cav- 
ity is operated in the TMno mode. The resonant frequency is 
2.8566 GHz and Q L = 14,700 + 50. The cavity is filled with 
gas (90% N2 and 10% iso-CztHio) with a gas-handling system. 
ISO-C4H10 is used as the quenching gas to remove background 
2 y-ray annihilation. 

3.1.5. Monte Cairo simulation 

The Monte Cairo simulation to use in analysis is performed 
using Geant4 lfl2ll . The low energy physics package PENE- 
LOPE IU3I1 is used and the geometry of the experimental setup 
is carefully input. The simulation is produced at each magnetic 
field strength and gas density. 

3.1.6. Data acquisition 

The prototype run was performed from 2 July 2009 to 24 
September 2009 using the large bore magnet with no compensa- 
tion (compensation magnets to reduce the uniformity to O(ppm) 
are planned but are not yet installed). In the overall period, the 
trigger rate was about 3.6 kHz and the data acquisition rate was 
about 650 Hz. The data acquisition was performed using NIM 
and CAMAC system. 

The trigger signal is the coincidence signal of yS-tagging sys- 
tem and y-ray detectors. Timing information of all PMT are ob- 
tained with a 2 GHz direct clock counting Time-to-Digital Con- 
verter (TDC: GNC-060) flU. A charge ADC (CAEN CI 205) 
is used to measure the energy information of the LaBr3 crystals 
while another charge ADC is used to measure the energy in- 
formation of the plastic scintillators. A crate controller (TOYO 
CC/NET) is used and obtained data are stored in a HDD of a 
Linux PC via Ethernet. The Zeeman transition has been mea- 
sured at various magnetic field strengths with a fixed RF fre- 



quency and power. The transition resonance lines are obtained 
at two gas densities (1.350 1 amagat and 0.891 6 amagat). 

3.2. Analysis 
3.2.1. Data analysis 

Figure [5] shows examples of measured timing spectra. The 
prompt peak coming from annihilation and |-) decay is fol- 
lowed by the decay curve of |+) and m z = ±1 states, then the 
constant accidental spectrum. A timing window of 35-155 ns 
is applied to select transition events. Figure [6] shows energy 
spectra, which are obtained by subtracting the accidental con- 
tribution using the timing window t =700-850 ns. 

The energy spectra are fitted using Monte Cairo simulation. 
The energy spectrum of RF-OFF 5 off can be expressed as 

5oFF=AoFF(5f y C + y 65f y C ) , (10) 

where Aoff is a normalizing constant, B is a 2y/3y decay ratio 
of RF-OFF, S is a Monte Cairo simulated energy spectrum 
of 2y decay, and 5 , ^ c is that of 3y decay. On the other hand, 
the energy specrum of RF-ON S on can be expressed as 

Son = Aon {s™ m + rS^) + A OFF (s™ at + BS^J , (H) 

where Aon is a normalizing constant and F is a 2y/3y decay 
ratio of RF-ON. The subscript "in" and "out" means the spec- 
trum of y-rays from Ps decays in the volume where RF power 
is applied or not applied, respectively. 

Typical fitted results are also shown in Figure [6] Figure [6(a)| 
is RF-OFF, and |6(b)| is RF-ON. The energy spectrum is mea- 
sured at different magnetic field strengths. The fitting range is 
370-545 keV. All the spectra are fitted successfully using MI- 
NUIT|Q. 

We take T—B as an amount of Zeeman transition. An example 
of resonance line obtained is shown in Figure Q Resonance 
lines can be fitted by S 2r /S 3r (RF-ON) - S 2r /S 3r (RF-OFF) • 
The theoretical function is calculated from Equations (0 and 
© numerically using RKF45 formula, with contributions from 
pick-off and slow positron. The typical Bo is 14.2 G. The RF- 
OFF function is obtained by substitute for Bq. Then the free 
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(a) Typical fitting of energy spectrum (RF-OFF). The black points 
are the data, the red line is 5^, c , the blue line is S^, c , and the pink 
line is sum of the Monte Cairo spectra. 
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(b) Typical fitting of energy spectrum (RF-ON). The black points 



are the data, the red line is 5^ lc : , 
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the aqua line is out , and the pink line is sum of 
the Monte Cairo spectra. Aoff and f} are fixed with the fitted values 
of RF-OFF, so that the free parameters are Aon and T. 



Figure 6: (Color online) Typical energy spectra and its fitting with Monte Cairo 
simulation. The data at 1.350 1 amagat, 0.865 816 9 T are shown and the spectra 
are normalized by the live time. 



RESONANCE LINE at 1.3501 amagat 




Table 1: Fitting result of the resonance lines. These uncertainties include errors 
from statistics of data, statistics of Monte Cairo simulation, uncertainty of RF 
power, and uncertainty of Ql value of the RF cavity. 



Gas density 


Ahfs 


Relative error 


^/ndf 


(amagat) 


(GHz) 


(ppm) 




1.3501(71) 


203.368 3(55) 


27 


0.910 


0.8916(23) 


203.3793(70) 


34 


0.483 



0.858 0.86 0.862 0.864 0.866 0.868 0.87 0.872 0.874 

MAGNETIC FIELD (T) 



parameters of fitting are Ahfs and C y . The fitting results are 
summarized in Table [TJ We use -33 ppm/amagat [1] as the 
material effect and obtain Ahfs in vacuum. 

3.2.2. Systematic Errors 

Systematic errors of the prototype run are summarized in Ta- 
ble |5] 

1 . Magnetic Field. The largest uncertainty in the prototype 
run is non-uniformity of the magnetic field. The weighted 
non-uniformity is 10.4 ppm. The offset and reproducibility 
of the magnetic field is measured to be 2 ppm. The calibra- 
tion uncertainty of NMR magnetometer is 1 ppm. These 
uncertainties are doubled because Ahfs is approximately 
proportional to square of the magnetic field strength. 

2. Monte Cairo Simulation. The magnetic field dependence 
of /3 is not exactly reproduced by the energy spectra fit- 
ting with Monte Cairo simulation. From the difference be- 
tween theoretical dependence and fitted result, this effect 
is estimated to be 18 ppm. The statistical uncertainty of 
Monte Cairo simulation is 17 ppm. 

3. RF System. The uncertainty of Qi value of RF cavity 
comes from reproducibility of the RF environment (~ 5%) 
and uncertainty of the Qi measurement method (0.6%). 
Its effect on Ahfs is estimated to be 6 ppm. RF power un- 
certainty comes from reproducibility of the RF environ- 
ment (~ 2%), uncertainty of power meters (0.05%), and 
uncertainty of the power meter calibration of temperature 
dependence (0.03%/K). It contributes to the error on Ahfs 
by 5 ppm. The uncertainty of RF frequency is 5 ppm. It 
directly affects Ahfs because Ahfs is approximately pro- 
portional to inverse of RF frequency. 

4. Material Effect. Thermalization of Ps can affect Ahfs by 
up to 20 ppm, but it has been not yet measured. Gas 
density dependence has not been measured in the proto- 
type run, so we have used the value from previous ex- 
periment. The uncertainty of the density dependence is 
4 ppm/amagat [1], and the uncertainty from iso-C 4 Hio is 
estimated to be less than 7.7 ppm/amagat. These result in 
7 ppm uncertainty of Ahfs ■ 

Other uncertainties are considered to be negligible. The sys- 
tematic errors discussed above are summed in quadrature. 



Figure 7: Resonance line at 1.350 1 amagat gas density. The circles and error 
bars are the data, and the solid line is the best fit result. The error bars include 
errors from statistics of data, statistics of Monte Cairo simulation, uncertainty 
of RF power, and uncertainty of Ql value of the RF cavity. 



3.3. Result 

The value of Ahfs obtained from the prototype run is 
Ahfs = 203.3804 ± 0.002 2(stat., 11 ppm) 
±0.008 l(sys.,40ppm)GHz, 



(12) 
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Table 2: Summary of Systematic errors of the prototype run. 
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7 
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which is consistent with both of the previous experimental val- 
ues and with the theoretical value. The uncertainty is an order 
of magnitude larger than the goal precision of O(ppm). 

3.4. Next steps 

The following improvements are planned for future measure- 
ments: 

1 . Compensation magnets will be installed and O(ppm) mag- 
netic field uniformity is expected to be achieved. 

2. The Monte Cairo simulation will be studied to reproduce 
the magnetic field dependence of energy spectra of Ps de- 
cays. Statistics of simulation will be reduced to O(ppm). 

3. The errors from RF system will be reduced to O(ppm) by 
carefully controlling the environment (especially the tem- 
perature) of the experiment. 

4. Measurements at various pressures of gas will be per- 
formed to estimate the material effect (the Stark effect). 
The accumulation of these measurements will result in an 
O(ppm) statistical error within a few years. 

5. The timing information allows for a measurement of Ps 
thermalization as a function of time 

HUH. 

We can thus 

precisely measure the material effect including the ther- 
malization effect. 



a few years which will be an independent check of the discrep- 
ancy between the present experimental values and the QED pre- 
diction. 
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4. Conclusion 

A new experiment to measure the Ps-HFS which re- 
duces possible common uncertainties in previous experiments 
has been constructed and the prototype run has been fin- 
ished. A value of A HF s = 203.3804 ± 0.002 2(stat.) ± 
0.008 l(sys.) GHz (41 ppm) has been obtained, which is consis- 
tent with both of the previous experimental values and with the 
theoretical calculation. Development of compensation magnets 
is underway with a view to obtaining O(ppm) magnetic field 
homogeneity for the final run. The final run will start soon. A 
new result with an accuracy of O(ppm) will be obtained within 
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